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Abstract 

Diatonns and other chlorophyll-c containing, or chromalveolate, algae are among the most productive and diverse 
phytoplankton in the ocean. Evolution arily, chlorophyll-c algae are linked through common, although not necessarily 
monophyletic, acquisition of plastid endosymbionts of red as well as most likely green algal origin. There is also strong 
evidence for a relatively high level of lineage-specific bacterial gene acquisition within chromalveolates. Therefore, analyses of 
gene content and derivation in chromalveolate taxa have indicated particularly diverse origins of their overall gene repertoire. 
As a single group of functionally related enzymes spanning two distinct gene families, fructose 1,6-bisphosphate aldolases 
(FBAs) illustrate the influence on core biochemical pathways of specific evolutionary associations among diatoms and other 
chromalveolates with various plastid-bearing and bacterial endosymbionts. Protein localization and activity, gene expression, 
and phylogenetic analyses indicate that the pennate diatom Phaeodactylum tricomutum contains five FBA genes with very 
little overall functional overlap. Three P. tricomutum FBAs, one class 1 and two class 11, are plastid localized, and each appears 
to have a distinct evolutionary origin as well as function. Class 1 plastid FBA appears to have been acquired by chromalveolates 
from a red algal endosymbiont, whereas one copy of class 11 plastid FBA is likely to have originated from an ancient green algal 
endosymbiont. The other copy appears to be the result of a chromalveolate-specific gene duplication. Plastid FBA 1 and 
chromalveolate-specific class 11 plastid FBA are localized in the pyrenoid region of the chloroplast where they are associated 
with (3-carbonic anhydrase, which is known to play a significant role in regulation of the diatom carbon concentrating 
mechanism. The two pyrenoid-associated FBAs are distinguished by contrasting gene expression profiles under nutrient 
limiting compared with optimal CO2 fixation conditions, suggestive of a distinct specialized function for each. Cytosolically 
localized FBAs in P. tricomutum likely play a role in glycolysis and cytoskeleton function and seem to have originated from the 
stramenopile host cell and from diatom-specific bacterial gene transfer, respectively. 

Key words: carbon concentrating mechanism (CCM), carbon metabolism, carbonic anhydrase, diatom, fructose 
bisphosphate aldolase, pyrenoid. 



Introduction 

It is believed that the diatom lineage emerged less than 300 
Ma. Considering the diversity of extant diatom species, 
comparable to that of angiosperms (Round et al. 1990), di- 
atoms have displayed a remarkable diversification rate and 
plasticity for adaptation to new environments. Mecha- 
nisms for assimilation and recycling of major nutrients, 
such as carbon (C), nitrogen (N), phosphorus (P), or iron 
(Fe) likely present notable differences from those known in 
green algae and plants and may provide insights about the 
ability of diatoms to thrive and diversify in aquatic environ- 
ments (Wilhelm et al. 2006; Armbrust 2009; Bowler et al. 



2010) . For example, diatoms possess a urea cycle, previously 
known only in metazoans, through which organic nitrogen 
may be recycled (Armbrust et al. 2004; Allen et al. 2006, 

2011) , and they may utilize siderophore-based iron uptake 
like in cyanobacteria (Allen et al. 2008). 

Diatoms are responsible for up to 40% of annual primary 
productivity in the ocean (Field et al. 1998; Granum et al. 
2005). Diatom RuBisCO (ribulose 1,5'bisphosphate carboxyl- 
ase/oxygenase) enzymes, however, have half saturation con- 
stants for CO2 of 30-60 |iM (Badger et al. 1998), whereas 
seawater typically contains around 10 |iM CO2, which implies 
that marine diatoms would be C02'limited in the absence of 
a carbon concentrating mechanism (CCM) (Riebesell et al. 
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1993) (Hopkinson et al. 2011; Reinfelder 2011). The relatively 
high affinity of diatom cells for CO2, coupled with the obser- 
vation that diatom populations commonly reach bloom den- 
sities, indicates that diatoms possess efficient CCMs (Badger 
et al. 1998; Roberts et al. 2007a, 2007b; Bowler et al. 2010; 
Hopkinson et al. 2011; Reinfelder 2011). 

In recent years, the CCMs, through which diatoms en- 
sure the supply of CO2 to RuBisCO, have attracted in- 
creasing attention (Roberts et al. 2007b). Green algae 
typically utilize bicarbonate transporters in combination 
with plastid, cytosol, and cell surface carbonic anhydrases 
(CAs), which catalyze conversion between CO2 and bicar- 
bonate, as a mechanism of promoting CO2 supply to Ru- 
BisCO, much like the biophysical CCMs that have been 
reported in cyanobacteria (Kaplan and Reinhold 1999). 
Diatom genome sequences unambiguously indicate the 
presence of different HCO3- transport systems (Kroth 
et al. 2008). In Phaeodactylum tricomutum, at least one 
such transporter has an N-terminal plastid-targeting pre- 
cursor (Allen et al. 2008; Kroth et al. 2008). Phaeodactylum 
tricorutum is also known to contain a plastid-targeted 
(3-carbonic anhydrase ((3'CA) (Montsant et al. 2005; Kitao 
et al. 2008; Tachibana et al. 2011). 

In diatoms, a putative C4-like mechanism has also been 
proposed to explain their efficient C fixation under CO2 and 
Zn limitation (Reinfelder et al. 2000, 2004; Reinfelder 2011). 
However, conflicting experimental data for support of C4 
photosynthesis in diatoms has been reported (Roberts 
et al. 2007a; McGinn and Morel 2008), and genomic data 
does not fully clarify the presence and localization of the en- 
zymes that may drive this mechanism (Kroth et al. 2008; 
Parker et al. 2008; Bowler et al. 2010). It is possible that in- 
dividual diatom species may rely to a different extent on 
both biophysical and biochemical CCMs in order to opti- 
mally regulate efficient inorganic carbon acquisition. 

Potential links between diatom CCMs and other aspects 
of C metabolism such as the direction of C flow between C3 
and Ce sugar pools and sustained regeneration of the C5 
RuBisCO substrate ribulose 1,5-bisphosphate (Ru 1,5-BP) 
have not been addressed. |3-CA and putatively plastidic 
fructose 1,6-bisphosphate aldolase (FBA) are downregu- 
lated and upregulated, respectively, in iron-limited P. tricor- 
nutum cells (Allen et al. 2008), which are known to display 
compromised photosystem reaction centers, reduced pho- 
tosynthetic electron transfer rates, decreased reductant pro- 
duction, and an inability to efficiently process absorbed 
photons (Behrenfeld et al. 1996, 2006; Milligan and Harrison 
2000). This has prompted speculation that plastidic (3-CA 
and FBA are part of a coordinated effort to regulate carbon 
flux within the plastid during periods of energy limitation 
(Allen et al. 2008). FBA genes have been noted for their pe- 
culiar phylogenetic distribution and highlighted in several 
comparative genomic studies (Kroth et al. 2005, 2008; 
Montsant et al. 2005) and due to important roles in balanc- 
ing Cg and C3 sugar pools in the Calvin-Benson cycle and 
glycolysis, as likely key regulators for the flow of small organic 
C molecules and CO2 in the plastid and cytosol. There are 
two evolutionarily unrelated FBAs, termed class I and class II, 



with a very complex phylogenetic distribution. These two 
types do not share sequence similarity with each other 
and their catalytic mechanism is different; unlike class I FBAs, 
class II FBAs are dependent on divalent cations and, there- 
fore, constitute a case of convergent functional evolution 
(Marsh and Lebherz 1992). 

Class I and class II FBAs catalyze the interconversion be- 
tween C3/C4 molecules and Cs/Cy molecules (Flechner 
et al. 1999). Therefore, FBAs can be involved both in carbon 
fixation and in glycolysis. In Calvin-Benson cycle reactions, 
the C5 molecule ribulose 1,5-bisphosphate (Ru 1,5-BP) is 
carboxylated and split into C3 molecules that are reduced 
to glyceraldehyde 3-phosphate (G3P). The formation of the 
hexose sugars glucose and fructose as well as regeneration 
of Ru 1,5-BP requires the action of FBA, which catalyzes the 
formation of the compound fructose 1,6-bisphosphate 
from two C3 G3P molecules and dihydroxyacetone phos- 
phate (DHAP). On the other hand, in glycolytic reactions 
FBAs split Cs molecules into C3 molecules, which are trans- 
formed into pyruvate. Pyruvate can then be transformed to 
acetyl coenzyme A and oxidized to CO2 through the Krebs 
cycle in the mitochondria or enter biosynthetic pathways 
such as fatty acid or amino acid biosynthesis. 

In animals, only class I FBAs are known, whereas fungi ap- 
pear to rely solely on FBA II (Jacobshagen and Schnarrenberger 
1990; Marsh and Lebherz 1992; Pelzer-reith et al. 1993) (table 
1). Red algae and glaucocystophytes have been reported to 
have class II FBAs in their cytosol and class I FBA in their 
plastids, although the sampling of these eukaryotic lineages 
has been limited to date (Anita 1967; Gross et al. 1994, 1999). 
The two classes of FBA have been detected in Eubacteria 
(including Cyanobacteria), although these organisms most 
typically utilize class II FBAs (Sanchez et al. 2002). Although 
FBA I activity was detected in the diatom P. tricomutum in 
an early study (Anita 1967), secondary endosymbiotic algae, 
commonly referred to as Chromalveolata (Cavalier-Smith 
2000), such as diatoms, dinoflagellates, haptophytes, and 
cryptophytes, which are believed to have inherited their plas- 
tids from red algae in a unifying single event (Yoon et al. 
2002), are generally thought to utilize class II FBAs both 
in cytosol and in plastids (Anita 1967; Rogers and Keeling 
2004). The universal occurrence of class II type A FBA in chro- 
malveolate plastids coupled with an apparent lack of red 
alga-like plastid-targeted FBA I was previously interpreted 
as evidence for a chromalveolate-specific gene replacement 
that supports the single origin of chromalveolate plastids 
(Patron et al. 2004). More recently, statistical analyses of 
phylogenomic data have been used to argue for falsification 
of the single origin of Chromalveolata in favor of more 
complex evolutionary scenarios such as serial acquisition 
of secondary plastids by distantly related hosts (Baurain 
et al. 2010). Whether chromalveolate plastids evolved from 
a single endosymbiotic event or from serial transfer be- 
tween diverse hosts, genomic data for key outgroups indi- 
cates an unambiguously common line of descent for 
Chromalveolata plastids (Janouskovec et al. 2010). 

Although a single gene family such as Fba will not resolve 
the conflicts surrounding the origin of chromalveolate plastids. 
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Table 1. Distribution of FBA Genes in Eukaryotes. 



Eukaryotic Lineage 


Eukaryotic Crown Group 


Plastid FBA 


Cytosol FBA 


References 


Plants 


Viridiplantae 


1 


1, IIB 


Jacobshagen and Schnarrenberger 










(1990); Pelzer-Reith et al. (1993) 


Red algae 


Viridiplantae 


1 


l^ 11 


Anita (1967); Kroth et al. (2005); This work 


Glaucophytes 


Viridiplantae 


IIB 


l^ 11 


Gross et al. (1994); This work 


Prasinophytes 


Viridiplantae 


l^ 11^ 


l^ (+ \^)^ 


This work 


Oomycetes 


Heterokonts 




II 


Kroth et al. (2005) 


Diatoms*^ 


Heterokonts 


(|po^ II, (iiPT 


11, l^ (+ 1'') 


Patron et al. (2004); Kroth et al. 










(2005); Montsant et al. (2005) 


Haptophyta 


Chromalveolates'* 




l^ 11^ 


This work 


Cryptophytes 


Chromalveolates'* 


1, (i|PT 




Patron et al. (2004); This work 


Apicomplexa 


Alveolates 




1 


Rogers and Keeling (2004) 


Dinoflagellates 


Alveolates 


(|Py)^ (llPyja 


l^ 11^ 


Kroth et al. (2005); This work 


Chlorarachniophyte 


Rhizaria 


1 


1 11^ 


Rogers and Keeling (2004) 


Animals 


Opisthokonts 




1 


Marsh and Lebherz (1992) 


Fungi 


Opisthokonts 




II 


Marsh and Lebherz (1992) 



Note. — The two known types of FBAs can participate in glycolytic reactions in the cytosol and in Calvin-Benson cycle reactions in the plastid. The two types are also widespread 
between Archaea and Eubacteria. indicates plastidic FBAll or FBA! localized to the pyrenoid. Pyrenoid localization of FBA has only been confirmed in P. tricornutum. 
^ Denotes first report in this work. 

^ This cytosolic type-! FBA is a divergent, bacterial-derived variant dissimilar to type-1 FBAs known in other eukaryotes. 

^ Phaeodactylum tricornutum has three putative plastid FBAs (the respective genes termed FBAC1, FBAC2, and FBAC5) and two cytosolic FBAs (FBA3 and FBA4). FBA4 and 
FBAC5 are type-1 FBAs. 

^ Haptophya and Cryptophyta are closely related to Heterokonts and together with diatoms and brown algae (stramenopiles), and dinoflagellates (alveolates) are known as 
chromalveolates. 



infornnation related to the phylogenetic distribution and 
function of different FBA types will help to clarify the evo- 
lutionary and physiological significance of different FBAs 
within and among different algal groups. Preliminary com- 
parative analyses of the distribution of FBA family genes in 
diatoms based on the genome sequence of the centric 
diatom Thalassiosira pseudonana and expressed sequence 
tag (EST) data from P. tricornutum revealed several surpris- 
ing features (Montsant et al. 2005). On the one hand, dia- 
toms appeared to be the only eukaryotes to possess 
a typically bacterial class I FBA, termed FBA4 (Kroth et al. 
2005), for which no targeting sequence was detected (i.e., 
likely a cytosolic enzyme). Additionally, the pennate dia- 
tom P. tricornutum appeared to have a plastid-localized 
class I FBA, FBAC5 (Kroth et al. 2005; Montsant et al. 
2005), which are usually found in the green and red algal 
lineages rather than in chromalveolate algae (Patron et al. 

2004) . This phylogenetic puzzle poses functional questions 
when considering that both T. pseudonana and P. tricornu- 
tum also have the class II plastid and cytosolic FBAs ex- 
pected in chromist algae (Patron et al. 2004), denoted 
FBAC1, FBAC2, and FBA3 in P. tricornutum (Kroth et al. 

2005) . We have therefore characterized each P. tricornutum 
FBA phylogenetically, using genome sequence or EST data 
now available for each of the major lineages of eukaryotic 
photosynthetic algae. To generate further information 
about functionality, we have also examined their subcel- 
lular localization and evaluated gene expression under 
different conditions. We specifically evaluated FBA tran- 
scriptional responses to Fe limitation because previous 
studies based primarily on EST sequencing suggested that 
certain FBA genes are responsive to Fe availability (Allen 
et al. 2006; Maheswari et al. 2010). Iron availability is also 
known to heavily influence diatom physiology and me- 
tabolism (Kustka et al. 2007; Allen et al. 2008) and distri- 
bution (Behrenfeld et al. 2006). 



Materials and Methods 

Phaeodactylum tricornutum Cultures 
Axenic cultures of P. tricornutum Bohlin clone CCMP2560 
(Provasoli-Guillard National Center for Culture of Marine Phy- 
toplankton, Bigelow) were grown in f/2 medium made with 
0.2 lim filtered and autoclaved seawater supplemented with 
f/2 vitamins and inorganic nutrients (Guillard 1975). Cell cul- 
tures were grown to midexponential phase at a temperature 
of 18 °C and a 16:8 light (150 |iE/m/sec):dark regime. Sterility 
was monitored by inoculation into peptone-enriched media 
to detect bacterial growth. For iron limitation experiments, Fe 
was precomplexed with EDTA (1:1.1 mobmol) and added 
(5 nM Fe) to Fe-free media. Cells were cultured in semicon- 
tinuous mode until growth rates stabilized. Cultures were di- 
luted with fresh medium at a rate of approximately 0.1-0.2 
per day depending on growth rate and were maintained in 
a state of exponential growth. Fe', the sum of all unchelated 
species, calculated according to equations of Sunda and 
Huntsman (2003) was (Fe' = 0.0024 Fej) under our condi- 
tions. Therefore, the Fe-limiting condition had an Fe' level of 
12.0 pmol liters"^ Fe' and the Fe replete condition had an 
Fe' value of 2.6 nmol-liters~^ Fe'. 

The availability of free Fe ions in the culture media was 
effectively buffered to maintain consistent Fe' levels over 
the course of the experiments. Growth rates were moni- 
tored by cell counting via microscopy before and after di- 
lution. Generally, the growth rate of cells in iron-limited 
cultures was around 60% of the maximum growth rate 
in iron-replete cultures. Iron limitation was confirmed by 
following the recovery of high growth rate following addi- 
tion of iron to a small volume of the culture. 

Cloning of P. tricornutum Genes 
The five Fba genes and a plastid-localized (3-CA (CA7) de- 
scribed earlier (Tanaka et al. 2005) were cloned into Gateway 
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diatom expression vectors (Siaut et al. 2007). The primers 
and cDNA templates used for polymerase chain reaction 
(PCR) amplification are indicated in supplementary table 
SI (Supplementary Material online). Correctly sized ampli- 
cons were cloned into Gateway pENTR/D-TOPO vectors (In- 
vitrogen) and sequenced. All FBA pENTR clones were then 
recombined into a diatom adapted pDEST_5' enhanced yel- 
low fluorescent protein (EYFP) through LR recombination 
reactions (LR clonase mix, Invitrogen), to generate pEXPR_ 
gene_EYFP vectors driving expression of the genes of interest 
with an EYFP tag at their Otermini (Siaut et al. 2007). FBA4, 
the bacteria-like class 1 FBA, was cloned both with and with- 
out a stop codon and recombined into pDEST vectors to 
generate N- and C-terminal EYFP fusions. The P-CA1 was 
used as a plastid marker for colocalization with FBA proteins. 
This gene was PCR amplified from the plasmid pre138- 
mptcal -enhanced green fluorescent protein (kindly pro- 
vided by Prof Yusuke Matsuda, Kwansei-Gakuin University, 
Hyogo, Japan), cloned into a pENTR vector, and recombined 
into C-terminal cyan fluorescent protein (CFP) fusion pDEST 
vectors to generate pEXPR_PtCA1_CFP. 

Transformation of Wild-Type P. tricornutum Cells 
The pEXPR plasmids (3 |ig) described above were each mixed 
with pAF6 plasmid (3 |ig), which drives expression of the 
phleomycin-resistance ShBle gene and precipitated onto 
tungsten microparticles (M-17, Biorad) of 1.1 jim in diameter. 
The constructs were transfected into P. tricornutum wild-type 
cells by helium-accelerated microparticle bombardment us- 
ing a BioRad PDS 1000 device following standard procedures 
described previously (Falciatore et al. 1999). Transformants 
were selected on solid medium (f/2 with agarose) containing 
100 |ig/ml phleomycin during 2 weeks, and individual colo- 
nies were then transferred to liquid medium for screening for 
EYFP signal. EYFP positive transformants were maintained in 
liquid media in the presence of the phleomycin analog xeocin 
(50 |ig/ml); which is lethal for wild-type cells. The EYFP pos- 
itive transformants have not lost their YFP signal over a period 
of 3-5 years, and the cultures have maintained resistance in- 
dicating long-term stability of transformed constructs. 

Epifluorescence Microscopy 

Transformant clones were screened for fluorescence 2 days 
after transfer to antibiotic-containing liquid medium. A 
Zeiss Axioskop equipped with an HBO 50W UV lamp 
was used, and the filter 93/XF104-2 (OMEGA), with exci- 
tation at 500 nm and emission at 545 nm, was used to de- 
tect YFP signal. The filter 93/XF102-2 (excitation 560 nm, 
emission 595 nm, OMEGA) was also used to acquire chlo- 
rophyll autofluorescence images. 

Confocal Microscopy 

EYFP-positive P. tricornutum clones detected by epifluores- 
cence microscopy were observed using a Leica TCS SP2 
confocal laser scanning microscope using an HCX PL 
APO IOOx/1. 40-0.70 oil objective. EYFP and chlorophyll 
were excited at 514 nm, CFP was excited at 458 nm, 
and the resulting fluorescence was filtered using a beam 



splitter (DD458/514). Fluorescence was detected using 
a photomultiplier tube with a bandwidth of 525-571 
nm for EYFP, 465-551 nm for CFP, and 620-710 nm for 
chlorophyll autofluorescence. Images were captured with 
the Leica confocal software and processed with the ImageJ 
1.34 software (Wayne Rasband, NIH). 

Transmission Electron Microscopy 
Three hundred microliters of an FBAC5-EYFP clone batch cul- 
ture were grown to late exponential phase and harvested by 
gentle centrifugation. Cells were fixed in 2% glutaraldehyde 
and embedded in LRWhite medium grade resin (London 
Resin Company). Ultrathin sections were treated with rabbit 
anti-GFP antiserum (AbCam) followed by goat anti-rabbit 
conjugated to electron-dense 10-nm gold particles as de- 
scribed previously (Lichtle and McKay 1992). Sections were 
observed with a Jeol CX2 electron microscope at 80 kV. 

Quantitative RT-PCR 

Cells were cultured in duplicate under Fe replete and deplete 
conditions, and mRNA was extracted from the four pellets 
(TRIzol, Invitrogen). RNA concentration was measured by 
means of an ND-1000 Spectrophotometer (Nanodrop), and 
Superscript 111 reverse polymerase (Invitrogen) was used to 
convert 100-ng total RNA into cDNA by reverse transcriptase- 
polymerase chain reaction (RT-PCR). The Fba gene tran- 
scripts were quantified by qPCR along with 18S rRNA as 
an internal control, in triplicate, with each reaction (25 
|il) containing 1 jil of the cDNA preparation, 200 nM of for- 
ward and reverse primers (supplementary table S2, Supple- 
mentary Material online), and 2x SYBR Green 1 PCR Master 
Mix (Applied Biosystems, CA). Reactions were run in an Op- 
ticon Chromo4 Mj Research Thermal Cycler (Bio-Rad, 
Hercules, CA). The cycling conditions comprised 10-min 
polymerase activation at 95 °C and 40 cycles at 95 °C for 
15 s and 60 °C for 60 s. The reaction was ended with 
a 5-min elongation step at 72 °C. Amplicon dissociation 
(melting) curves were recorded after cycle 40 by heating from 
60 to 95 °C with a ramp speed of 0.5 °C every second. The 
Fba threshold cycles (Cj) were normalized to three endog- 
enous control gene 18S rRNA, TATA binding protein, and H4 
(Siaut et al. 2007). Fold difference between Fe replenished 
and Fe starved cultures measured by means of the 2~^^ 
^T method (Livak and Schmittgen 2001), where AACj 

~ (Cj, target ~ ^T, 185 )5 nM ~ (Cj, target ~ ^T, 185 )l |aM- 

FBA Activity Assays 

Protein concentrations of cellular extracts were determined 
by Bradford assay (BioRad, Hercules) according to protocols 
of the manufacturer. Enzyme activity was measured in an 
optimized coupled test based on the protocol by Gross 
et al. (1999). The standard reaction assay contained 100 
mM Tricine, pH 7.5; 10 mM MgCl2; 1 U/ml triosephosphate 
isomerase (TIM); 1 U/ml glycerol-3-phosphate dehydro- 
genase; NADH (0.2 mM); cellular extract (varying vol- 
umes). The decreasing concentration of NADH was 
monitored spectrophotometrically at 2 = 366 nm. The 
reaction was started at 25 °C by the addition of 2 mM 
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fructose'1,6'bisphosphate. Cellular extracts were preincu- 
bated in the presence of bivalent ions (10 mM either AAgCl2, 
MnCl2, ZnCl2) or 10 mM EDTA on ice for 1 h. Class II FBAs 
are dependent on divalent metal ions and are, therefore, in- 
hibited by EDTA (Perham 1990). 

Phylogenetic Inference 

For each of the diatom FBA genes, the predicted protein 
model was searched for homologs in a comprehensive ge- 
nomic database, which we assembled from the nonredun- 
dant GenBank CDS database (nr), JGI microbial genomes, 
and EST libraries (NCBI; ESTdb) for lineages that do not 
have complete genome data (e.g., red algae and dinoflagel- 
lates). Using the iTree phylogenomic pipeline (Moustafa 
et al. 2010), the homologs for the queries were consolidated 
at each of the taxonomic ranks (i.e., species, genus, class, 
order, phylum, and kingdom) to maximize taxon coverage, 
while maintaining a practical number of taxa. Each query 
and its homologs were aligned using MAFFT (Katoh et al. 
2002; Armbrust et al. 2004; Moustafa et al. 2010) with up to 
1,000 refinement iterations. Multiple sequence alignments 
were manually examined and refined to ensure that key 
taxa were not missing and to remove ambiguous sites from 
the alignments. Next, molecular phylogenies were inferred 
using PhyML (Guindon et al. 2009) with a maximum likeli- 
hood and LG amino acid substitution model (Le and Gascuel 
2008) and discrete gamma distribution (Yang 1994). 
Branch confidence values were estimated using the ap- 
proximate likelihood-ratio test (Anisimova and Gascuel 
2006) implementation in PhyML. Only support values 
greater than 50% were indicated on the graphical repre- 
sentation of the phylogenetic trees. Accession numbers 
for all sequences analyzed are given in supplementary 
table S3 (Supplementary Material online). 

Results and Discussion 

Intraceiiuiar Locaiization and Phytogeny of P. 
tricomutum FBAs 
Plastidic FBAs 

Observation of multiple positive transformants of C- 
terminally tagged FBAC1, FBAC2, and FBAC5 confirmed that 
all three of these proteins are localized in the plastid. The 
FBAC1-EYFP fusion appears as small spots (usually only 
one or two) inside the plastid, in a very discrete and precise 
fashion (fig. 1 A), with virtually no EYFP background observed 
elsewhere in the plastid or the cytosol. Although, phyloge- 
netically, FBAC2 seems to be a duplication of FBAC1 (Kroth 
et al. 2005), its localization was not like FBAC1 but rather 
appeared to be diffuse throughout the plastid (fig. IB). 
On the other hand, the class I plastid FBA, FBAC5, typical 
of higher plants and algae with primary plastids, formed 
well-defined spots similar to FBAC1-YFP (fig. 1C). 

To further investigate the subplastidal localization of 
FBAC1 and FBAC5, P. tricomutum P-type CA, PtCAl, which 
was recently reported to localize to the pyrenoid (Tachibana 
et al. 2011) was cloned and tagged with a C-terminal CFP 
fusion, and cotransformed with YFP-tagged FBAC1 and 




Fig. 1. Localization of plastid FBA-EYFP fusion proteins. (A) FBAC1- 
EYFP, (B) FBAC2-EYFP, (C) FBAC5-EYFP. On the left light microscopy 
image, chlorophyll autofluorescence is shown in red and YFP 
fluorescence in green. Bars, 2 jam. 

FBAC5. The CA1-CFP signal indeed appeared as punctate 
spots within the plastid, which colocalized with both 
FBAC1 and FBAC5 (fig. 2A and B). It therefore appears 
that all three proteins localize to the same intraplastidial 
compartment. 

Immunolocalization of EYFP in transformants expressing 
FBAC5-EYFP showed clearly that the protein localizes to 
the pyrenoid (fig. 3). We can conclude that PtCAl, FBAC1, 
and FBAC5 all localize to pyrenoids within diatom plastids. 
Pyrenoids are electron-dense bodies known to contain 
a protein matrix that consists mainly of RuBisCO (Griffiths 
1970). Although some classes of microalgae, including cer- 
tain dinoflagellates, have been reported to activate a CCM 
in the absence of an observable pyrenoid (Morita et al. 
1998; Ratti et al. 2007), the pyrenoid is generally considered 
to be an important factor for CCM function when present; 
however, precise function remains unclear (Giordano et al. 
2005). Notably, enlargement of pyrenoids has been ob- 
served under CO2 limitation (Izumo et al. 2007), and model 
simulations of mass spectrometric measurements of cellu- 
lar carbon fluxes in P. tricomutum wild-type and P-CA over 
expression cells suggested that the pyrenoid is the likely site 
of CO2 elevation to concentrations required to saturate 
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Fig. 2. Coiocaiization of p last id FBA-YFP with CA1-CFP. (A) FBACl- 
EYFP + CA1-CFP, (B) FBAC5-EYFP + CA1-CFP. Chlorophyll 
autofluorescence is shown in red, YFP fluorescence in green, and CFP 
fluorescence in blue. Left image: CFP fluorescence, center: YFP 
fluorescence, right image: combined fluorescence. Bars, 2 |im. 



carbon fixation (Hopkinson et al. 2011). Although the bio- 
chemical significance of FBA proteins in the diatom pyre- 
noid remains to be investigated, it is clear from these 
experiments that certain members of plastidic class 1 
and 11 FBA in diatoms have undergone functional diversi- 
fication related to partitioning some aspects of Calvin cycle 
activity into the pyrenoid; near the site of CO2 fixation. 

Phylogenetic topology for plastid-localized FBAC5 (fig. 4) 
suggests a monophyletic relationship among chromalveolate 
plastids and indicates inheritance from the red algal endo- 
symbiotic ancestor. It is interesting to note that, compared 
with previous analyses of class 1 FBA distribution, more thor- 
ough taxonomic sampling suggests that cytosolic FBAC5 or- 
thologs are highly diversified and present in most major 
eukaryotic groups including all the chromalveolates as well 
as Rhizaria and glaucophytes. Photosynthetic strameonpiles 
such as the multicellular brown algae Ectocarpus siliceous 
and the pennate diatom Fragilariopsis cyimdrus also ap- 
pear to contain cytosolic class 1 FBA. Curiously, one new 
well-supported cytosolic clade includes cryptophytes. 



dinoflagellates, and haptophytes and appears to have origi- 
nated in these taxa via endosymbiotic gene transfer (EGT). 

The plastidic FBAC5 clade is also much expanded rela- 
tive to previous reports. The most likely explanation for 
plastidic FBA 1 (FBAC5) occurrence in the diatom, hapto- 
phyte, dinoflagellate algae clade is gene duplication within 
Archaeplastida prior to the diverge of red algae, followed by 
EGT from red algae to chromalveolates coupled to selective 
gene loss within some heterokonts (e.g., centric diatoms). 
Haptophytes and Ectocarpus appear to have a second copy 
of plastidic FBAC5 that resides in a clade that includes Rhi- 
zaria and freshwater green algae. 

Class 11 A FBAs include the plastid-localized copies FBAC1 
(pyrenoid) and C2 and cytosolic FBA3. FBAC1 and C2 are 
members of highly resolved clades comprised of all the major 
marine phytoplankton taxa including prasinophytes, heter- 
okonts, dinoflagellates, haptophytes, and cryptophytes 
(fig. 5). The occurrence of prasinophytes within the FBAC2 
clade is a new and significant finding and provides some 
insights into the mystery surrounding the origin of these 
genes in chromalveolates (Kroth et al. 2005). Considering 
the well-supported affiliation between prasinophytes and 
chromalveolates, it seems unlikely that chromalveolate 
plastid-localized FBAC2 arose by duplication of class 11 cyto- 
solic FBA3. Rather a more plausible scenario is that, like class 
1 plastid-localized FBAC5, class 11 plastid FBAC2 was also ac- 
quired via EGT. Plastid-targeted nuclear genes of green algal 
derivation within chromalveolates have been proposed to 
have arisen from a cryptic endosymbiont related to a 
prasinophyte-like green alga (Moustafa et al. 2009). An 
alternative interpretation is that such genes have simply 
been lost in extant red algae. Regardless of whether the 
source of the FBA 11 plastidial EGT event was a green 
or red algae, it seems apparent that the gene was acquired 
from a member of the Archaeplastida lineage. 

The phylogenetic relationship between prasinophytes 
and chromalveolates suggests that pyrenoid-associated 
FBAC1 is the product of an ancient chromalveolate-specific 
duplication of FbaC2. Due to the occurrence of FBAC1 in 
each of the four major chromalveolate lineages (hetero- 
konts, alveolates, cryptophytes, and haptophytes), this 




Fig. 3. Localization of Phaeodactylum tricomutum FBAC5 by transmission electron microscopy. A clone expressing the FBAC5_EYFP fusions 
was fixed, and the intracellular EYFP was immunolocalized by means of an anti-GFP antibody. The antibody bound strongly to an intraplastidial 
structure in FBAC5_EYFP transformant cells. Chloroplasts (C) and pyrenoids (Pyr) and gold particles detected with electron microscopy are 
indicated. Unlabeled cells did not show any labeling (data not shown). 
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Fig. 4. Protein maximum likelihood phylogeny of class 1 FBA. Numbers at the nodes indicate approximate likelihood-ratio test support values. 



duplication nnust have occurred prior to chronnalveolate 
diversification. FBAC2, on the other hand, appears to be 
present only in haptophytes and diatoms and is likely to 
have been lost within at least some chromalveolates. 

Cytosolic FBAs 

Cytosolic class II FBA, termed FBA3, displayed a typically 
cytosolic localization in YFP fusion experiments (fig. 6). Part 
of the YFP fluorescence appears to concentrate around the 
plasma membrane or as thin strands within the cell; these 
regions of dense signal could represent cytosol-rich areas 
between vacuolar compartments and other organelles 
(as observed in clones expressing YFP without any targeting 
signal; data not shown). 

The bacterial-class I FBA, called FBA4, showed a similar 
cytosolic localization both for N-terminal (data not shown) 
and C'terminal EYFP fusions (fig. 7). In some cells, quite 
finely defined intracellular lines were observed, much more 
clearly than in the case of the FBA3. The conspicuous and 
precisely defined strands surrounding the plastid in the 
FBA4 transformants may suggests interactions with the cy- 
toskeleton. Interaction of FBA with the cytoskeleton may 
relate to a role of FBA in motility. In the alveolate Plasmo- 
dium falciparum, cytosolic FBAs have actin binding activity 



and play a role in coordinating interactions between cell 
surface adhesion and cytoskeleton control of parasite mo- 
tility (Jewett and Sibley 2003). 

Prior to diatom genome sequencing, FBA4 had only been 
described in bacteria. All three diatom genomes sequenced 
to date contain Fba4, and it was also found in ESTs from the 
cold-adapted polar diatom Chaetoceros neogracile (Jung et al. 
2007) (fig. 8). The pelagophyte Aureococcus and prasino- 
phyte Micromonas also encode genes for FBA4 that form 
a well-supported clade distinct from diatoms (fig. 8). 

Expression and Activity of P. tricornutum FBAs 
A set of 15 cDNA libraries obtained from mRNA derived 
from cell cultures grown in different conditions were gener- 
ated in the context of the P. tricornutum genome sequencing 
project (Bowler et al. 2008; Maheswari et al. 2010). Between 
4,000 and 10,000 ESTs were sequenced from each of these 
cDNA libraries, giving a total collection of around 132,000 
sequences that were assembled into about 12,000 transcrip- 
tional units (Maheswari et al. 2010). See Maheswari et al. 
(2010) for a description of culture conditions. 

All five FBAs are expressed in several libraries (fig. 9), ex- 
cept for plastid class I FBA, FBAC5, which appears to be 
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transcribed only under iron limitation (FL). Similar to 
FBAC5, plastid class II FBA, FBAC1, was also upregulated 
in Fe limited cultures and additionally in response to blue 
light (BL). The other plastid class II FBA, FBAC2, was also 
abundant in the blue light library but not under iron lim- 
itation. The blue light library was constructed from material 
harvested after exposure to 2 h of blue light following 48 h 
in the dark and therefore should represent conditions fa- 
vorable for C fixation, protein synthesis, and rapid growth. 



Cytosolic FBAs also showed an interesting differential ex- 
pression profile. The class II FBA FBA3 showed clear upre- 
gulation in blue light and to a lesser extent in NO3' and Fe 
limitation. The bacterial-type FBA, FBA4, was expressed un- 
der a range of different conditions. 

The observed EST patterns under Fe limitation were 
confirmed by quantitative RT-PCR (table 2). It was appar- 
ent that plastid class I FBA, FBACS, and plastid class II FBA, 
FBAC1, are upregulated up to 90- and 35-fold, respectively. 




Fig. 6. Localization of FBA3-EYFP fusion proteins. On the left light microscopical image center: YFP fluorescence in yellow and chlorophyll 
autofluorescence is shown in red. YFP fluorescence is shown in green on the right in order to more easily distinguish the chloroplast. Bars, 2 jam. 
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Fig. 7. Localization of FBA4. (A) FBA4-EYFP fusion; YFP fluorescence is shown in yellow. (B) EYFP-FBA4 fusion; YFP fluorescence is shown in 
green, chlorophyll autofluorescence in red. Bars, 2 lam. 



whereas the other plastid class II FBA, ^^kCl is strongly 
downregulated under Fe limitation. The cytosolic class II 
FBA, ^^k?>, was also found to be more highly expressed un- 
der iron limitation (over lO-fold difference), whereas the 
expression of the bacterial class I V^M is unresponsive 
to changes in Fe concentration. Furthermore, in experi- 
ments using P. tricornutum strain Pt4 (UTEX646), we found 
a clear increase of total FBA and class I FBA activity (mea- 
sured after inhibition of class II aldolases) compared with 
cells grown under Fe replete conditions (fig. 10). 

Conclusions 

FBAs are central to many aspects of cellular C metabolism 
and are likely relevant for diatom ecophysiological mecha- 
nisms of managing acclimation to nutrient limitation and 
other factors, such as temperature and light availability, that 
influence C fixation and turnover. Diatoms have a multiplic- 
ity of dynamically expressed plastid and cytosolic FBA en- 
zymes likely adapted for management of specific anabolic 
and catabolic functions in different intracellular compart- 
ments under shifting environmental conditions. Preliminary 
investigation of FBA function in the pennate diatom P. tri- 
cornutum suggests little functional overlap. The two types of 
cytosolic FBAs, for example, have very different expression 
profiles. Diatoms have a class II cytosolic FBA (FBA3), which 
is a divalent cation-dependent enzyme, as well as a bacterial- 
like cytosolic class I variant (FBA4) unknown in eukaryotes 
other than diatoms. The pattern of expression for FBAS, en- 
coding a class II cytosolic FBA, suggests a role in both syn- 
thesis and mobilization of storage polysaccharides in 
response to the onset of light and nutrient (Fe and N) lim- 
itation conditions, respectively. On the other hand, the ex- 



pression profile of the bacterial- like FBA4 indicates 
a relatively high level of expression in the triradiate morpho- 
type and sensitivity to the presence of silicic acid. Silicate 
availability greatly influences the distribution and abundance 
of diatoms in nature, and upregulation in triradiate morpho- 
type cells may be suggestive of a role in cytoskeletal dynam- 
ics. Experimental data indicate that FBA4 expression does 
not change in response to Fe availability, which combined 
with the EST data does not suggest a role in central carbon 
metabolism (e.g., Calvin-Benson cycle and glycolysis). The 
well-defined cytosolic strands suggestive of cytoskeletal in- 
teraction are not as readily observed in FBA3-EYFP. The pre- 
cise function of diatom-specific FBA4 remains to be 
uncovered; however, localization and expression data are 
suggestive of a diatom-specific functional specialization 
rather than a role in conventional C metabolism. 

Whereas the pyrenoid-associated plastid-localized class I 
FBA, FBAC5, may be upregulated and active specifically un- 
der the low Calvin-Benson cycle C flux conditions of Fe lim- 
itation, class II pyrenoid-associated FBA (FBAC1) was 
upregulated under Fe limitation as well as in response to 
the onset of blue light, following a 48-h dark period. By con- 
trast, non pyrenoid-associated plastid class II FBA {FBAC2) 
was downregulated in response to Fe limitation and stimu- 
lated by blue light. Upregulation of FBAO and FBAC2 gene 
expression in response to the onset of light suggests a role for 
these genes under favorable conditions for C fixation. 

Confirmation of pyrenoid localization for (3-CA in P. tri- 
cornutum (Tachibana et al. 201 1 ) supports the idea of a bio- 
physical CCM based on transport of HCO3' into the plastid 
followed by conversion of HCO3- from the stroma to CO2 
in the pyrenoid near the active site of carbon fixation 
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(Hopkinson et al. 2011), as previously hypothesized 
(Roberts et al. 2007a, 2007b). The occurrence of class I 
and II FBAs, FBAC5 and FBAC1, in the pyrenoid of P. tricor- 
nutum might implicate linkage between some components 
of the Calvin-Benson cycle and the CCAA. FBAs catalyze the 
first step in the regeneration of the 5 carbon (C) ribulose 
compound that is the substrate for RuBisCO, CO2 fixation, 
and generation of new C3 sugars. It is possible, therefore, that 
the FBA catalyzed conversion of C3 to molecules serves to 
partially regulate CO2 fluxes into the Calvin-Benson cycle. In 
the vascular plant Arabidopsis thaliana, FBA is 1 of 3 Calvin 
cycle enzymes, along with RuBisCO and SBPAse, that was 
most sensitive to biological perturbations and found to have 
an important rate-limiting role in regulating the carbon 
assimilation flux of the Calvin cycle (Haake et al. 1998; Sun 
et al. 2003). 

In photosynthetically active plastids, FBA activity is nec- 
essary for regeneration of Ru 1,5-BP and should therefore 
operate principally in the direction of Fru1,6-BP formation, 
as suggested by antisense inhibition of FBA in potato plants 
(Haake et al. 1998). Under optimal conditions for C fixation, 
prior to generation of C5 RuBisCO substrate, com- 
pounds, generated from FBA activity, prime the regenera- 



tion phase of the Calvin-Benson cycle, and subsequently 
for every molecule that is used, 5 of 6 are recycled 
and 1 of 6 constitutes net gain. However, in situations 
in which inorganic C fixation is constrained, the reverse re- 
action (Zgiby et al. 2000) may occur to permit sufficient 
flow of C through the triose pool for respiration and bio- 
synthesis of other C molecules. Indeed, gene expression, 
metabolite, and physiological data indicate an overall re- 
trenchment of C fixation in response to Fe limitation, a con- 
dition which has been globally assayed in detail in P. 
tricornutum (Allen et al. 2008), including downregulation 
for genes encoding P-CA, phosphoribulokinase (PRK), 
and fructose 1,6-bisphosphatase (FBPase) (Allen et al. 
2008), all of which would be coherent with reduction of 
C5 regeneration. Consistent with the idea of FBA-mediated 
production and associated chloroplast export of trioses, 
plastid triose isomerases and plastid-envelope triose/hex- 
ose phosphate translocators are upregulated under Fe lim- 
itation (Allen et al. 2008). 

In some microalgae, storage polysaccharides have been 
shown to be mobilized in nitrogen (N), C, Fe, and light- 
limited cultures (Varum and Myklestad 1984; Granum and 
Myklestad 2001; Allen et al. 2008), and in addition to 
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Fig. 9. Expression of Phaeodactylum thcomutum FBA genes according 
to their normalized frequency of observation in a range of cDNA 
libraries representing different environmental conditions (Maheswari 
et al. 2010). Normalized expression values were obtained by dividing 
the number of FBA ESTs in the different libraries by the total number 
of sequences generated for the different libraries. The R value was 
defined in Stekel et al. (2000), an R value above 12 indicates 
nonspurious differences in the EST frequency across libraries (i.e., 
differential expression). The features of the 15 libraries have been 
published elsewhere (Allen et al. 2008; Sapriel et al. 2009; Maheswari 
et al. 2010). Si- = f/2 with artificial seawater, no Si added; Si+ = 350 
|iM metasilicate; Oval = Oval strain CCAP1052/1B; Triradiate = 
Triradiate strain NEPCC640; Blue light = Cells grown in the dark 
for 48 h and subjected to blue light for 1 h; Tropical = Tropical strain 
CCMP633 grown at a suboptimal temperature of 18 °C; N replete = 
1.12 mM chemostat culture; N03 starved = 50 |iM in chemostat 
culture for 12 days; Urea adapted = 50 jiM Urea; Fe deplete = 5 nM 
Iron; Low DD = 0.5 |ig/ml 2E,4E-decadienal for 6 h; High DD = 5 iig/ 
ml 2E,4E-decadienal for 6 h (for a full description of culture 
conditions, see Maheswari et al. 2010). 

supplying respiratory substrate, may be imported into the 
plastid to ensure the continued biosynthesis of essential 
metabolites (such as fatty acids) under impaired C fixation 
conditions. Consequently, a possible upregulation of genes 
encoding FBAC1 and FBAC5 could imply a role in the main- 
tenance of carbon flow through C3 and C4 pools for the 
synthesis of lipids, nucleotides, and other metabolites. 
Given that diatoms do not have an oxidative pentose phos- 
phate pathway in their plastids (Michels et al. 2005) and 
neither do they have storage polysaccharides within them, 
the mobilization of C3 or molecules across plastid mem- 
branes during impaired C fixation could be necessary for 
the continued supply of C building blocks for anabolic pro- 
Table 2. Fold Variation of FBA Gene Expression Observed in Fe- 
Limited Cultures with Respect to Fe Replete Cultures (5 nM vs. 10 
liM). 





Experiment 1 


Experiment 2 


FbaCl 


26.1 ± 2.2 


45.3 ± 5.2 


FbaC2 


-22.9 ± 0.7 


-36.0 ± 33.5 


Fba3 


13.0 ± 3.1 


16.7 ± 2.9 


Fba4 


1.2 ± 0.2 


1.4 ± 0.2 


FbaC5 


86.4 ± 23.9 


100.8 ± 28.0 



Note. — Two identical experiments were conducted on different dates; the fold 
change was calculated by the 2~^^ method (see Materials and Methods). The 
values shown are averages and standard deviations of two replicate mRNA 
extractions. 
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Fig. 10. FBA activity in Phaeodactylum thcomutum (UTEX646) cells 
grown for 5 days in Fe deplete and replete medium. Enzymatic 
activity is given as mU/mg total cellular soluble protein. The mea- 
surements were performed in either presence or absence of EDTA to 
measure class I FBA activity or total FBA activity, respectively. 

cesses. In contrast, the expression profile of cytosolic class II 
FBA, FBA 3, which is upregulated under N and Fe limitation 
and in response to the onset of blue light, is consistent with 
a role in catabolic C flow through glycolysis. 

Most algal genomes encode between 2 and 3 plastid lo- 
calized FBAs and 1 or 2 cytosolic copies. In diatoms, class II 
cytosolic FBA, FBA3, appears to be the only FBA that is 
likely to have originated with the secondary endosymbiotic 
host nucleus; the other four diatom FBAs appear to have 
originated from lateral transfer from bacteria and from 
plastid acquisition. In contrast to plants and freshwater 
algae, all the major marine unicellular eukaryotic phyto- 
plankton lineages, prasinophytes, haptophytes, dinoflagel- 
lates, diatoms, and cryptophytes, encode plastid-targeted 
copies of class I and II FBA. Chromalveolate taxa are distin- 
guished by a second copy of class II plastid FBA, FBAC1, 
which, in P. tricornutum, appears to function in the pyre- 
noid region of the plastid. The most parsimonious expla- 
nation for the origin of FBAC1 is a chromalveolate-specific 
duplication of FBAC2, which appears to have originated by 
EGT from green algae. The overall pattern of FBA diversity 
within and between chromalveolates could be interpreted 
with the context of either the chromalveolate or serial hy- 
pothesis for the origin of plastids in chromalveolates (Baurain 
et al. 2010). In the latter scenario, FBAC2 duplication would 
have occurred in the founder alga that subsequently engaged 
in eukaryote-eukaryote endosymbiosis with various hosts. 
The presence of FBAC2 in diatoms and haptophytes perhaps 
suggests shared host ancestry at the time of FBA2 acquisi- 
tion — via EGT from a cryptic green algal endosymbiont. 
The idea that chromalveolates have recruited genes sepa- 
rately from red and green algae (Moustafa et al. 2009) seems 
to be typified by plastid-localized FBAs whereby one copy, 
class I FBA (FBAC5), appears to have originated from red algal 
EGT and the other two copies, class II FBAs (FBAC1 and 
FBAC2), seem to have been acquired via EGT from a prasino- 
phyte-like green alga followed by a chromalveolate-specific 
duplication. 
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Supplementary tables SI -S3 are available at Molecular Biology 
and Evolution online (http://www.mbe.oxfordjournals.org/). 
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